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This  addendum  to  Draper  Laboratory  Report  R-977,  Inertial 
Navigation  System  Standardized  Software  Development,  is  the 
first  in  a series  of  informal  memos  documenting  work  completed 
since  May,  1976.  For  convenience  it  has  been  published  and 
bound  with  a cover. 
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SHORT  TERM  FOLLOW-ON 
TO  STANDARD  SOFTWARE  PROGRAM 


(MAY  - JULY  1976/ 


Background 

By  the  end  of  April  1976,  the  Numerical  Simulator,  NUMSIM 
program  had  been  checked  out  and  was  running  (with  AFAL's 
Profile  Generator  program,  PROFGEN)  on  the  AFAL  CDC  6600. 

The  final  report  had  been  rough  drafted  and  was  submitted  to 
documentation. 

However,  the  need  for  some  modifications  to  PROFGEN  and 
hence  to  NUMSIM,  were  evident.  In  addition,  some  post-process- 
ing programs  were  needed  for  efficient  use  of  the  PROFGEN/NUMSIM 
package,  and  a short  parametric  study,  and  some  long  runs  with 
representative  hardware/software  configurations  were  desired. 

The  short  term  follow-on  efforts  were  directed  towards 
fulfilling  these  objectives  (as  specified  in  a Statement  of 
Work  and  Task  Breakdown) . 

Summary 

1)  AFAL  modified  PROFGEN  to  compute  and  output  the  inte- 
grals of  specific  force,  in  each  of  the  local  vertical, 
space  stable,  and  strapdown  frames  - rather  than  the 
specific  forces  per  se  as  had  been  provided  previously. 

2)  CSDL  modified  NUMSIM  to  accept  these  integrals  rather 
than  instantaneous  values  of  specific  force  with  subse- 
quent (approximate)  integration  as  had  previously  been 
the  case. 

3)  NUMSIM  was  modified  to  compute  the  sum  of  the 
squares  of  the  errors  over  its  output  cycle  time.  This 
provided  the  capability  of  generating  short  or  long  term 
rms  errors  in  post-processing  programs  (in  addition  to 
instantaneous  errors) . 
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4)  A plotting  program,  to  provide  the  capability  of 
generating  CALCOMP  plots  of  any  of  the  errors  in  (3) , 
was  added. 

5)  Checkout  runs  employing  the  modified  PROFGEN  and 
NUMSIM  programs  indicated  that  the  errors  in  NUMSIM  out- 
puts for  gimballed  systems  (local  vertical  or  space 
stable)  could  be  made  arbitrarily  small.  Small  errors 
in  strapdown  system  outputs  occur  principally  due  to 
discontinuities  (steps)  in  the  angular  rates  between 
PROFGEN  outputs. 

6)  The  variable  precision  version  of  NUMSIM  VUMSIM  was 
checked  out  on  the  CDC  6600. 

7)  A short  parametric  study  was  performed  to  obtain  pre- 
liminary estimates  of  the  sensitivities  of  system  per- 
formance to  changes  in  hardware/software  configuration. 

8)  Long  runs  were  performed  using  the  same  mission 
profile  and  hardware  configuration  but  varying  only  the 
software.  Generally,  the  simpler  (or  "baseline")  soft- 
ware was  shown  to  be  adequate  for  a moderate  accuracy, 
aircraft  INS.  An  anomaly  in  the  performance  of  the 
"upgraded"  software  during  benign  flight  conditions 
led  to  continued  investigations  during  July. 

9)  The  investigation  revealed  that  a 24-bit  mantissa 
was  not  adequate  for  the  position  computations  (direc- 
tion cosine  matrix,  update,  orthonormalization)  in  an 
accurate  INS.  The  mantissa  was  extended  for  these  compu- 
tations-to-simulate-32-bit,  fixed  point  computations. 

This  eliminated  the  anomaly.  (A  similar  precision  exten- 
sion would  be  useful  in  the  velocity  summing  computations.) 
Sign  errors  (in  both  PROFGEN  and  NUMSIM)  on  one  component 
cf  gravity  were  corrected. 

10)  More  detailed  descriptions  of  the  short  term  follow- 
on  efforts  and  results  are  presented  in  the  body  of  this 
report. 
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I.  Changes  Made  to  Simulator 

A.  Used  revised  TAPE 20  input  and  input  format  as  follows: 

i)  Expanded  "TRAJIN"  COMMON  block  in  all  routines  to 
include  the  three  sets  of  specific  forces.  These  values  are 
stored  in  the  SFIS  array  as  indicated  on  'Revised  TAPE 20 
Format' . 

ii)  Block  data  routine  modified  to  initialize  SFIS  to 
zero. 

iii)  INREC  routine  modified  to  read  the  updated  (23 
word)  records. 

iv)  CMINTG  was  modified  so  that  when  IPC(24)  is  not 
equal  to  zero,  the  delta  velocities  are  found  by  summing 
the  PROFGEN  specific  force  integrals,  as  picked  up  from 
SF$T . 

v)  Changes  were  made  to  SSINTG,  LLINTG,  and  SDINTG  so 
that  when  IPC(24)  is  not  equal  to  zero,  each  routine 

a)  picks  up  information  about  the  integral  of  spe- 
cific force  from  SFIS, 

b)  '■.orrects  it  to  NUMSIM  frame,  misaligning  it  (if 
appropriate)  for  the  local  level  case,  and 

c)  stores  it  in  SF$T  for  use  by  CMINTG. 

B.  Generated  statistics  and  wrote  them  (using  revised 
TAPE 30  output  format)  and  printed  them  as  follows: 

i)  Executive  was  modified  so  that  at  the  end  of  each 
navigation  cycle 

a)  OUTUNI  was  called, 

b)  If  IPC(22)  is  not  equal  to  zero,  a n ew  subroutine 
SSQALL  is  called. 

ii)  A new  COMMON  block  SUMSQS  was  added  to  Block  Data, 
PLTAPE , PRINTR  and  SSQALL.  The  block  contains  SUMSQl(lO) 
and  SUMSQ2 (10) . SUMSQ1  contains  the  sums  of  the  squares 
of  the  difference, (between  PROFGEN  and  NUMSIM  at  the  end 
of  navigation  cycles) , since  the  beginning  of  the  NUMSIM 
run.  The  differences  are  of  the  following  ten  items  in 
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this  order: 
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1. 

Latitude 

- degrees 

2. 

Longitude 

- degrees 

3. 

Alpha 

- degrees 

4. 

Altitude 

- feet 

5. 

Velocity  up 

- feet/sec 

6 . 

Velocity  east 

- feet/sec 

7. 

Velocity  north 

- feet/sec 

8. 

Roll 

- degrees 

9. 

Pitch 

- degrees 

10. 

Heading 

- degrees 

SUMSQ2  has  the  sum  of  the  squares  as  above,  only  accumulated 
since  the  last  tape  output. 

iii)  Output  COMMON  block  was  redone  so  that  DLAT,  DLONG , 
DALF,  DH,  DV ( 3 ) , DETA(3)  appear  in  the  order  shown  in  this 
sentence. 

iv)  PRINTR  was  expanded  so  that,  when  IPC(22)  is  not 
equal  to  zero,  the  square  roots  of  the  values  in  SUMSQ1 
are  printed  out. 

v)  PLTAPE  was  expanded  so  that  it  writes  out  informa- 
tion in  the  new  TAPE 30  format,  which  implies  taking  the 
mean  of  the  values  in  SUMSQ2,  writing  out  the  mean,  and 
then  clearing  SUMSQ2. 


I 

m 


1 


1 
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Miscellaneous  Cleanups,  Corrections  and  Improvements 

1)  In  the  exec,  PTIME  and  PLTIME  were  added  to  the 
SIMPAR  NAMELIST.  This  was  done  to  allow  easier  control 
over  starting  both  printout  and  tape  output. 

2)  PRINTER  and  PLTAPE  were  fixed  to  agree  with  the 
current  initialization  scheme.  PLTAPE  was  changed  to 
"rewind"  its  output  file  during  initialization. 

3)  A test  for  STH  equal  to  zero  was  added  to  SDATUD  to 
prevent  a possible  division  by  zero. 
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Reviseu  "TAPE 30"  Format 


Record  Length 
1 10 


Contents 

40  character  title,  stored  four 
characters  per  word 


2-n  21 

1) 

time,  in  seconds 

2) 

A 

latitude 

3) 

A 

longitude 

where  A is 

4) 

A 

alpha 

difference,  at 

5) 

A 

altitude 

' time' seconds , 

6) 

A 

velocity 

(up) 

between  PROFGEN 

7) 

A 

velocity 

(east) 

and  NUMSIM; 

8) 

A 

veloci ty 

(north) 

9) 

A 

roll. 

10) 

A 

pitch 

ID 

A 

heading 

and  where  Z is  the 

12) 

Z 

latitude 

mean  sum  of  the 

13) 

Z 

longitude 

squares  of  the 

in 

z 

aplha 

differences  (at 

15) 

z 

al titude 

the  end  of  each 

16) 

z 

velocity 

(up) 

navigation  cycle) 

17) 

z 

velocity 

(east) 

between  PROFGEN 

18) 

z 

velcoity 

(north) 

and  NUMSIM, 

19) 

z 

roll 

since  the  last 

20) 

z 

pitch 

tape  output. 

21) 

V 

o 

heading 
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Revised  "TAPE20"  Format  for  PROFGEN  Output 


1 


Record  # 
1 

2 

3 

4-n 


Contents 

Unmodified 

Unmodified 

Unmodified 

Words  1-14  are  unmodified  words, 

15-23  are  integrals  (over  the  output  period) 
of  specific  force  as  follows: 

15)  North  component  (for  a = 0)  x 
West 
Up 

Forward  i 
Right  Wing  I Strapdown 

Down 


16) 

17) 

18) 

19) 

20) 
21) 
22) 
23) 


Local 

Level 


North 

West 

Up 


Space  Stable 


* 


The  foi: owing  two  pages  define  the  PROFGEN  and  NUMSIM  coordinate 
conventions. 
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D.  Changes  to  NUMSId  to  implement  more  precision  in  DCMUPD 
and  ORTHO 

1362  WRITE  (6,691) 

1364  691  FORMAT  (*b  THIS  VERSION  HAS  GREATER  PRECISION", 
1366  IN  DCMUPD  AND  OP.THO"/"b  AND  ASSUMES  FRACTL  (=40') 

10092  COMMON (PRECIS) FRACTL,  EXPO 

10094  INTEGER  FRACTL,  EXPO 

10675  FRACTL  = FRACTL+8 

10745  FRACTL  = FRACTL- 8 

10755  FRACTL  = FRACTL+8 

10765  FRACTL  = FRACTL- 8 

14542  COMMON/PRECIS/FRACTL,  EXPO 

14544  INTEGER  FRACTL,  EXPO 

14555  FRACTL  = FRACTL+8 

14625  FRACTL  = FRACTL- 8 

E.  Fix  to  "North  component  of  gravity"  (actually  projection 

A 

of  this  onto  y axis  in  LLWA  frame) . In  subroutine  GRAV  change 

G$2P (2)  = -COEF*A0P$2P (3,2) 
to 

G$2P (2)  = COEF*A0P$2P ( 3, 2) 

11  Changes  Made  to  PROFGEN 

a)  By  AFAL;  calculations  of  integrals  of  specific  force 
(Note  that  IRITE  must  be  1 for  printed  output  to  be  correct) . 

b)  By  CSDL;  took  out 

CALL  ACCLRTN  (FX,  FY,  FZ) 

at  lines  11620  and  14400. 
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III . Summary  of  Short  Parametric  Study 

Fifteen  8-second  runs  were  performed  to  gain  some  insight 
into  the  sensitivity  of  the  navigation  errors  to: 

a)  word  length 

b)  navigation  computation  cycle 

c)  algorithm  "order" 

d)  quantization 

The  first  three  runs,  sequence  numbers  1,  2,  and  3,  used  the 
baseline  software  with  coarsest  quantization,  a 1/8  second  navi- 
gation computation  cycle,  and  the  full  precision  word  (corres- 
ponding to  48 -bit  mantissa)  and  represented  a strapdown  (SD) . 
a local  vertical  (LV) , and  a space  stable  (SS)  inertial  measuring 
unit  (IMU)  respectively. 

The  next  three  runs,  sequence  numbers  4,  5,  and  6,  differed 

from  the  first  three  only  in  the  reduction  of  the  mantissa  (frac- 

tional part)  of  the  word  length  to  24  bits. 

The  next  three  runs,  sequence  numbers  7,  8 and  9,  all  employ 
the  higher  order  algorithm,  with  the  full  word  length  and  no 
quantization,  for  a LV  IMU.  They  differ  only  in  the  frequency 
of  navigation  computations,  128,  32,  and  8 times  a second. 

The  next  three  runs,  sequence  numbers  10,  11,  and  12,  employ 

the  higher  order  algorithm,  with  full  wordlength,  8 Hz  navigation 
computations  for  a LV  xMU  with  coarse,  moderate  and  fine  quantiza- 
tion of  incremental  velocity  and  gimbal  angles.  The  values  of 
VQUANT  are  0.04,  0.01,  and  0.0025  feet  per  second  per  pulse,  while 
the  corresponding  values  of  AQUANT  are  0.384,  0.096,  and  0.024 
milliradians  respectively. 

The  final  three  runs,  sequence  numbers  13,  14,  and  15,  are  all 
8 Hz,  LV  IMU  cases  with  full  word  length  and  no  quantization.  They 
differ  only  in  the  portions  of  the  higher  order  algorithm  in- 
cluded, which  latter  are  controlled  by  the  flags,  IPC  (28),  IPC 
(29),  and  IPC  (30'.  With  all  3 flags  set  to  zero  (Seq.  #13)  the 
baseline  algorithm  is  mechanized.  With  IPCOO)  = 1,  (Seq.  #14) 
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the  exact  expression  for  angular  velocity  replaces  the  approxi- 
mate expression  used  in  the  baseline  case.  With  IPC(29)  = ipc  (30) 
= 1,  the  interpolation,  extrapolation  and  change  of  flow  of  the 
upgraded  algorithm  is  added  to  previous  case.  The  resultant  is 
the  higher  order  algorithm  with  the  exception  that  a first 
order  direction  cosine  matrix  update  is  performed. 

RESULTS  OF  RUNS,  SEQUENCE  NOS.  1 THROUGH  15 

In  most  cases,  the  output  of  each  run  was  displayed  on  the 
CRT  of  the  Hazeltine  2Q00  (72  columns) , and  printed  out  on  the 
80  column  satellite  roll  type  printer.  This  process  was  rather 
tedious  and  the  resulting  formats  were  ungainly  compared  with 
the  normal  line  printer  outputs  - so  these  printouts  are  not 
included. 

In  some  cases,  Calcomp  plots  were  produced  and  are  included. 

A complete  set  of  plots  consists  of  three  sets  of  ten  plots 
each  representing  the  instantaneous,  the  rms  values  (between  out- 
put intervals  and  the  time  rms  (from  start  of  run  to  indicated 
output  t me)  of  each  of  the  following  navigation  errors: 

- latitude,  longitude,  wander  angle  (in  degrees) 

- altitude  (in  feet) 

- up,  east,  and,  north  velocity  (in  feet  per  second) 

- roll,  pitch,  and  heading  ?in  degrees) 

The  more  significant  results  of  the  short  parametric  study 
are  summarized  in  Tables  1 through  6 in  declining  order  of  im- 
portance (for  the  particular  8 second  simulated  flight) . Note 
that  only  the  rms  errors  at  the  end  of  the  runs  are  compared. 

Calcomp  plots  of  the  instantaneous  errors  from  a strapdown 
run  preceding  Sequence  #1  are  shown  in  plots  #1  through  10.  This 
represents  best  available  strapdown  performance  (with  the  pre- 
sent PROFGEN  dynamics).  Hand  plots  of  the  east  and  north  veloc- 
ity errors,  every  1/128 th  of  a second,  are  included  to  illustrate 
the  error  build-up  during  the  l/2g  turn. 
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Calcomp  plots  of  the  short  term  runs  ("average”)  errors 
from  the  sequence  #1  "baseline"  strapdown  run  are  shown  in 
plots  #13  through  22. 

Hand  plots  of  the  rms  east  velocity  errors  during  the  se- 
quence #1,  2,  3 runs  show  the  velocity  errors  approaching  the 
asymptotes  expected  due  to  incremental  velocity  quantization 
(plot  #25). 
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Table  1 


Comparison  of  Time  RMS  Errors 
at  End  of  8 Sec,,  8 Hz  LV  Runs 
- Effects  of  Algorithm,  Quantization,  Word  Length 
Ideal  vs.  Baseline 


Error 

Type 

Magnitude 

Ratio 

Units 

Ideal 

Seq.#9 

Baseline/wQ 
Seq. #5 

#5/#  9 

Long. 

y deg. 

.00242 

8.827 

3640 

Lat. 

y deg. 

.0125 

16.46 

1320 

Alt. 

feet 

.00241 

.6384 

254 

V east 

mfps 

.0110 

19.13 

1740 

V north 

mfps 

.0138 

21.89 

1590 

V up 

mfps 

.0677 

21.47 

317 

Roll 

y deg. 

.0822 

4.034 

49,000 

Pitch 

y deg. 

.6734 

10.525 

15,700 

Heading 

y deg. 

.0017 

16.063 

9,450,000 
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Table  2 


Comparison  of  Time  RMS  Errors 


at  End  of  8 Sec, , 8 Hz,  LV  Runs 


Ideal  vs.  Baseline  Algorithms 
(No  quantization) 


Error 

Units 

Ideal 
' (Seq. #9) 

|-  - - 

♦Long . 

y deg 

.00242 

Lat. 

y deg 

.0125 

Alt. 

feet 

.00241 

V east 

m feet/sec 

.0110 

V north 

m feet/sec 

.0138 

V up 

m feet/sec 

.0677 

Roll 

y deg 

.0822 

Pitch 

y deg 

.6734 

Heading 

y deg 

.0017 

Baseline 
(Seq. #13) 

11.18 

.6310 

.6153 

.1152 

.1611 

2.946 

.6175 

8.438 

7.908 


Ratio 

#13/#9 


7.5 

12.5 

4650 


♦Velocity  error  during  1st  sec  of  east  acceleration  causes 
11.6  y deg  longitude  error. 


Table  3 


Error 

Type 


V east 

V north 

V up 


♦Roll 

Pitch 

Heading 


Comparison  of  Time  RMS  Errors 


at  End  of  8 Sec  . LV  Runs 


vs.  Navigation  Comp.  Cycle 


Nav.  Com 


Units 


y deg, 
y deg. 
feet 


32 

Seq. #8 


8 

Seq. #9 


.000007 

.000139 

.002424 

.000166 

.000742 

.012452 

.000010 

.000150 

.002407 

.001244 

.001390 

.011036 

.000062 

.000860 

.013759 

.000410 

.004324 

.067733 

.035517 

.036532 

.082230 

.493597 

.498660 

.674350 

.000033 

.000118 

.001736 

Ratios 

#8. 

/#7 

00 

19.9 

17.4 

4.5 

16.8 

15.0 

16.0 

♦Attitude  at  128  Hz  for  all  of  #7,  #8,  #9 
(278  ydeg.  = 1 sec) 

* 12.74  ydeg.  = 1 foot  (Lat.) 

I 3.88  ydeg.  = 1 foot  (Long,  at  45°  Lat.) 
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Error 

Type 


Table  5 

Comparison  of  Time  RMS  Errors 
at  End  of  8 Sec.,  8 Hz  LV  Runs 
Effects  of  Quantization  & All  Else 
vs.  Quantization  Only 


Magnitude 


Units 

U deg. 
U deg 
feet 


Seq.  #10 
.1582 
.2218 
.07443 


Sec.  #5 
8.827 
16.46 
.6384 


Ratio 

#5, 

'*10 

5.58 
74.3 

8.58 


V east 

V north 

V up 

mfps 

mfps 

mfps 

19.52 

22.56 

21.71 

19.13 

21.89 

21.47 

Roll 

u deg. 

2,736. 

.4,034. 

Pitch 

V deg. 

9,989. 

10,525. 

Heading 

U deg. 

8,068. 

16,063. 
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Lonq  Term  Simulations 


Background 


At  the  AFAL-CSDL  meeting  on  3th  and  9th  of  Jun  >t  was  agreed 
that  the  mission  to  be  used  for  the  long  runs  should  be  a three 
hour  F-4  combat  interdiction  mission  for  which  AFAL  had  the  PROFGEN 
inputs  on  file.  The  file  name  was  supplied  to  CSDL. 

The  characteristics  of  the  IMU  and  navigation  computer  were 
defined  by  AFAL  as  follows: 

1)  IMU  - local  wander  azimuth 

- 0.040  fps  per  pulse  incremental  velocity  quantization 

- 14-hit  gimbal  angle  encoders 

- 0.40  arc  seconds  per  pulse  gyro  torquer  quantization 

2)  Navigation  Computer 

-floating  point  (8  bit  exponent  and  24  bit  mantissa) 

- 0.25  second  navigation  (and  attitude)  iteration  rate. 

PROFGEN  was  to  be  run,  using  the  specified  mission,  at  16  Hz, 
and  the  version  which  computes  the  integrals  of  specific  force. 

Two  runs  were  to  be  performed:  one  using  the  "baseline"  soft- 
ware and  the  other  the  "upgraded"  or  "full"  software.  Error  out- 
puts would  be  printed  every  60  seconds  (including  the  rms  errors) . 
Calcomp  plots  would  be  made  including  the  difference  between  the 
instantaneous  errors. 

Results  would  be  provided  by  an  informal  memo  type  report 
(later  changed  to  an  addendum  to  the  Final  Report) . 
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Early  Results  and  Interpretation 

First  attempted  PROFGEN  run  aborted  after  800  seconds  of  simu- 
lated time.  It  was  also  noted  that  file  specified  by  AFAL  started 
not  at  the  beginning  of  the  mission,  but  at  the  start  of  the  dynamic 
phase  and  included  evasive  maneuvers,  attack,  evasive  maneuvers 
and  climb  out. 

Two  VUMSIh  runs  were  performed  using  the  "baseline"  and  the 
"upgraded"  software  (referred  to  as  sequence  numbers  16  and  jl7 
respectively) . The  instantaneous  and  time  rms  position  and  velo- 
city errors  were  plotted  at  60  second  intervals  from  the  Hazelrine 
outputs  and  are  shown  in  plots  26  through  37.  A full  set  of  30 
Calcomp  plots  for  the  "baseline"  software  run  (sequence  #16)  are 
included  as  plots  38  through  67. 

As  anticipated,  the  "upgraded"  software  case  provided 
smaller  position  errors  than  the  "baseline"  case. 

A rough  order  of  magnitude  estimate  of  the  relative  sizes  of 
the  errors  is  provided  by  scaling  of  the  error  plots  (see  Table  7) 
while  a more  precise  estimate  is  given  by  the  ratio  of  the  time 
rms  errors  (see  Table  8) . Attitude  errors  were  almost  identical 
in  either  case  and  reflected  the  expected  values  due  to  encoder 
quantization. 

The  most  significant  result  of  these  runs  was  the  growth  rc te 
of  the  horizontal  position  errors  in  the  "baseline"  case.  The  rate 
error  was  on  the  order  of  0.10  knots,  or  compatible  with  the  compu- 
tational portion  of  the  error  budget  for  a moderate  accuracy  air- 
craft inert-  1 navigation  system. 
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TABLE  7 


ROM  ESTIMATE 
of 

ERROR  RATIOS 
from 

PLOT  SCALING 


(Baseline/Upgradcd) 


Error 

Name 

Plot  Scale 
Ratio 

Latitude 

2.5 

Longitude 

2.5 

Attitude 

50.0 

North  Velocity 

1.0 

East  Velocity 

1.0 

Vertical  Velocity 

20.0 

TABLE  8 

UPGRADED/BASELINE  SOFTWARE 
Ratio  of  RMS  Errors  at  660  Sec 


Error 

Units 

Baseline 

Upgraded 

B/U 

Longitude 

U degrees 

241 

53 

4.55 

Latitutde 

U degrees 

215 

62 

3.47 

Altitude 

feet 

11.8 

0.22 

53.60 

East  Velocity 

nxfeet/scc 

42 

55 

0.765 

North  Velocity 

mfeet/sec 

30 

26 

1.15 

Up  Velocity 

mfeet/sec 

370 

15 

24.9 
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FIRST  LONGER  RUNS  - RESULTS  AND  INTERPRETATION 


A second  PR0FG3N  through  6300  + seconds  was  performed  (the  roll 
rate  used  in  horizontal  turns  was  reduced  from  450  to  250  degrees 
per  second . ) 

"Baseline"  and  "full"  software  VUMSIM  runs  were  performed 
(through  4800  seconds) . The  resulting  RMS  errors  at  4500  sec- 
onds are  summarized  in  Table-  9.  Calcomp  plots  of  the  difference 
of  the  instantaneous  errors  between  the  two  runs  are  shown  in 
plots  68  through  77. 

In  general,  the  navigation  errors  from  the  "full"  software 
case  were  much  smaller  than  the  corresponding  errors  from  the  "base- 
line" case.  The  anomaly  occurred  in  latitude  with  the  "full"  soft- 
ware run;  here  the  rim  latitude  error  was  about  four  times  larger 
than  in  the  "baseline"  case. 

The  appearance  of  this  anomaly  initiated  a month  long  investi- 
gation to  localize  and  identify  the  source  of  the  error  or  errors 
and  fix  same. 

Table  9 

RMS  Errors  (Seq.  #16,  17)  at  4500  Sec. 


Error 

Units 

BSW 

FSW 

Ratio  B/F 

Lat. 

y deg. 

390 

1525 

0.26 

Long. 

y deg. 

1127 

360 

3.13 

Alpha 

y deg. 

539 

262 

2.06 

Alt 

feet 

4.55 

0.09 

50.5 

Vx  (up) 

milli  fps 

144 

12.2 

11.8 

Vy  (east) 

milli  fps 

173 

44.9 

3.86 

V_  (north) 
z 

milli  fps 

131 

110 

1.19 

Roll 

milli  deg. 

5.7 

5.9 

0.97 

Pitch 

milli  deg. 

6.1 

6.1 

1.00 

Yaw 

milli  deg. 

12.3 

12.9 

0.95 
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Latitude  Anomaly  Investigation 


Examination  of  the  north  velocity  and  latitude  errors  indicated 
that  the  position  error  was  not  the  integral  of  the  velocity  error  in 
the  "full"  software  run  (sequence  #17) . 

Initial  Hypotheses 

Since  we  were  dealing  with  errors  at  this  time  it  was  not 
clear  whether  the  errors  were  due  to  some  inadequacy  in  the  variable 
precision  simulator  program  (VUMSIM)  or  in  the  profile,  i.e.,  incon- 
sistencies in  the  position,  velocity,  and  integrals  of  specific 
force  output  by  PROFGEN. 

Special  Tests  and  Results 


Numerous  special  simulation  runs  over  part  or  all  of  the 
trajectory  used  in  sequences  #16  and  #17,  were  performed  to  isolate 
the  problem.  Most  of  these  runs  employed  the  version  of  the  simu- 
lator without  variable  precision,  (NUMSIM)  with  no  quantization  and 
with  16  Hz  navigation  and  attitude  computations. 

All  these  runs  display  relatively  small  systematic  position 
and  velocity  errors,  the  position  errors  were  the  integrals  of  the 
velocity  errors  and  the  familiar  forms  of  the  Schuler  oscillations 
were  apparent,  e.g„,  the  velocity  errors  were  of  the  form  A sinu>st 
then  the  corresponding  position  error  was  of  the  form  (A/u>s) 

(1  - cosaj  t)  . The  magnitude  of  A corresponded  to  an  accelerometer 
bias  error  or  "tilt"  error  ranging  from  a fraction  of  a micro-g  or 
microradian  up  to  2 or  3 micro-g' s or  microradians. 

AFAL  acknowledged  the  existence  of  a sign  error  in  one  of 
the  horizontal  components  of  "g"  at  altitude,  which  would  account 
for  the  kind  and  magnitude  of  errors  mentioned  in  the  preceding 
paragraph. 
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Simulation  runs  switched  to  the  variable  precision  form 
(VUMSIM)  where  the  48-bit  case  corresponded  exactly  to  NUMSIM. 

The  32-bit  VUMSIM  results  departed  microscopically  from  the  48-bit 
case.  When  the  corresponding  24-bit  case  was  run,  the  anomaly 
obligingly  reappeared. 

Further  testing  and  analysis  indicated  that  24-bits  were 

too  few  for  position  integration.  This  situation  was  further 

aggravated  by  orthonomalization  which  turned  the  direction  cosine 

matrix  errors  into  rotational  errors  (which  latter  appeared  as 

- 3 

drift  rate  errors  on  the  order  of  8 x 10  degrees  per  hour) . 


Fixes  and  New  Results 

Since  32  bits  was  good  and  24  bits  was  bad,  VUMSIM  was  modi- 
fied to  perform  the  direction  cosine  matrix  update  and  orthonormali- 
zation using  8 bits  more  than  was  used  in  the  remainder  of  the  pro- 
gram. The  DCM  per  se  is  stored  with  the  added  8 bits. 

Runs  at  4 Hz  with  this  fix  did  not  display  the  large  latitude 
errors  characterizing  the  original  anomaly  (note  that  no  quantiza- 
tion was  included  at  this  point) — see  plots  #78  and  #79. 

Next,  the  sign  on  the  one  component  of  "g"  w as  changed  in 
PROFGEN  and  the  profile  was  regenerated.  A 16  Hz  NUMSIM  run  with 
the  "full"  software  and  no  quantization  revealed  that  the  longitude 
error  had  about  doubled  compared  with  the  "old"  PROFGEN  instead  of 
vanishing  (see  piots  #80  and  #81). 

After  performing  several  4 Hz  runs  with  and  without  quantiza- 
tion but  with  extended  precision  in  the  DCM  operations  (see  plots 
82,  83,  84,  85,  86  and  87)  with  rather  inconclusive  results,  the 
"gravity"  equations  for  NUMSIM/VUMSIM  were  reexamined  for  the  nth 
time.  The  suspect  level  component  of  "g"  was  found  to  have  been 
originally  compatible  with  PROFGEN  (both  had  the  wrong  sign 
initially)  so  this  sign  was  also  reversed  and  the  appropriate  4 Hz, 
24  bits  (32  for  DCM)  VUMSIM  runs  with  quantization  were  repeated 
with  results  not  significantly  different  from  those  shown  in  plots 
82,  83,  86  and  87. 
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This  concluded  the  computer  work  performed  to  date. 

PROFGEN  Outputs  Before  and  After  Change  in  Sign  of  GY 

Over  the  F-4  mission,  changing  the  sign  of  GY  in  PROFGEN 
introduces  change  in  east-west  acceleration  of  about  2 micro-g's 
over  most  of  the  flight.  Doubly  integrated,  over  5000  seconds, 
this  would  represent  about  3000  microdegiees  of  longitude  or 
about  800  feet.  Applied  to  a Schuler  loop,  the  peak  position  error 
(after  half  a Schuler  period)  would  be  about  1/10  of  the  above. 

The  actual  difference  in  the  PROFGEN  longitude  due  to  this 
change  in  acceleration,  at  4800  seconds  into  the  run,  was  1.2 
microdegrees  or  about  4 inches.  This  might  lead  one  to  conclude 
that  the  PROFGEN  position  (and  velocity)  are  not  quite  compatible 
with  its  specific  force  (integrals) . 

Conclusions  and  Recommendations 

1)  Extended  precision  is  required  for  the  position  integra- 
tion - 24  bits  is  not  adequate  while  32  bits  is  more  than  adequate. 
Therefore,  the  DCM  update  should  be  performed  using  32-bit  fixed 
point  arithmetic  (corresponding  to  the  present  VUMSIM  mechaniza- 
tion) . 

2)  Similar  extended  precision  would  be  useful  in  the  velocity 
summing  block.  (This  is  not  currently  mechanized  in  VUMSIM) . 

3)  The  use  of  PROFGEN  should  be  restricted  to  providing  an 
incremental  velocity  profile.  The  reference  data  should  be  supplied 
by  the  "best"  software  NUMSIM  run,  i.e.,  full  word  length,  no  quanti- 
sation, short  computation  cycle,  etc.  (The  existing  "differencing" 
program  will  accommodate  this.) 

4)  Effects  of  quantization  should  be  reevaluated  after  incor- 
poration of  extended  precision  velocity  summing. 

5)  The  effects  of  step  inputs  from  PROFGEN  becloud  the 
propagation  of  computational  errors  during  benign  flight  conditions. 
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6)  The  throughput  capabilities  of  the  existing  remote  ter-  , t 
minal  are  not  conducive  to  efficient  analysis  of  the  data  available  | 
at  the  central  processor.  1 
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